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Abstract. The purpose of these lectures is to give a pedagogical overview of cosmological pertur- 
bation theory, following the lectures given during the school. The topics treated are: 

I - The background 

II - Scalar/Vector/Tensor decomposition and the gauge issue 

III - The example of the tensor modes 

IV - Density fluctuations, transfer function and power spectrum 

V - Initial condition theory: quantum vacuum fluctuations 

Most of the material presented here is available in many well-written reviews or textbooks, so in 
order to avoid unnecessary heavy presentation as well as to make sure I forget nobody, I will only 
cite the review paper [1] as well as the book [2] from which most of the figures have been taken. 
Useful extra information and different perspectives can be also found in [3] (in particular the review 
articles by A. Linde on inflation, J. Martin on the quantum aspect of initial condition and their 
subsequent squeezed evolution and C. Ringeval on the numerical evolution of perturbations). Ref. 
[4] provides a personal vision of S. Weinberg with many original proofs to well-known results, 
and [5] describes in more details the relevant physics for calculating the quantities actually to be 
compared with the data. Finally, all numerical figures are taken from the Particle Data Group [6] 
whose latest update is always available on the linked site. 



Even though these lectures concern cosmological perturbation theory, I felt an introduc- 
tion to the background could be welcome, would it be only to fix the notations, set the 
framework and make apparent what the problems and questions are. 

Cosmology is the part of physics that studies the Universe as a whole, trying to 
make models of its overall evolution and its structure. As such, it is a quite peculiar 
branch of physics, as by definition there is only one Universe - hence the name - 
and it is impossible to make any experiment on either its evolution or structure! From 
these considerations, we immediately see that cosmology will be endowed with various 
intrinsic limitations which I will discuss in due turn. 

How do we, practically, describe cosmology? To begin with, one needs a theoretical 
framework providing the evolution equations. This will be general relativity: the Uni- 
verse will be seen as a 4-dimensional manifold, space-time, endowed with a metric g^ v 
whose dynamics follows from Einstein equations 
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INTRODUCTION: THE BACKGROUND 




(1) 



where c, the velocity of light in vacuum, will be set to unity in all further calculations 
(along with h where it should have appeared in the final section of these lectures), G N is 
Newton's constant, is the stress-energy tensor of the matter - discussed later - and 
A the cosmological constant. We know from observations that the latter is probably not 
vanishing, contrary to what was supposed until recently, but we can however consider 
its influence as another matter fluid and include it in T^ v . Therefore, one can, without 
lack of generality, send A —> in Eq. (1). 

The Einstein tensor = R^ v - \g^ v R is defined in terms of the Ricci tensor R^ Y and 
scalar R = R^ = g^ v R^ v , the former stemming from a contraction of the Riemann tensor 
through R^ v = R"^. Finally, the relation with the metric itself is made with the definition 

and the Christoffel symbols are given in terms of the metric by 

r % ~ ( d agv/3 + dpg va - d v g a fj) . (3) 

This completes the geometrical explanation, i.e. the left hand side of Eq. (1). The next 
question, more physical in a way, now is: what is the matter content of the Universe? 
With this content well-defined, one can in principle find the relevant solutions of Einstein 
equations. General solutions of these equations are of course not known, so a less 
ambitious program consists in trying to find out a simple model for which we do have 
solutions! For this, we will need to impose some constraints. 

Before we even embark into describing the model itself, let us mention that we need 
to confront the following limitations: 

- The Universe is unique by definition, so the usual methodology of physics is not 
applicable as we can neither compare with other similar objects to evaluate how generic 
our observations are nor redo experiments! 

- We are observing the Universe from a single location in both space and time that we 
did not choose. In particular, this implies a question about the history of the Universe 
and the specific moment we happen to observe it. 

- Observations, as it turns out, are limited to our backward light cone, see figure 1. 

- For a given set of data, there are possibly many space-times corresponding to the 
observations. Again, as we have only one set of observations and since we cannot redo 
the experiment consisting in having the Universe evolving from the Big-Bang to now, 
we have no way to make sure our interpretation of the data is the correct one. We need to 
make some hypothesis on the nature of the structure of space-time and verify those. Only 
the large number of repetitions of observational data can reduce the risk of confusion 
between different models. Nowadays, we have so many data explained by one single 
model that it has become extremely difficult, if not altogether impossible, to come up 
with a different but equally successful model. One can already notice, at this point, that 
a further complication with any model of the Universe is that most of the hypothesis are 
hard to verify, since we can actually model our Universe, i.e. the observable one. The 
description of the actual Universe, which may even be infinite in size, does not belong to 
the realm of physics as we will never have access to it. Indeed, the Universe is probably 
much larger than the observable Universe over which we can collect data. 
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FIGURE 1. Limitations due to our given and unchangeable position in the Universe. We can only 
measure objects located inside our past lightcone, would they be geophysical data on our worldline or 
astrophysical data seen though its light emission, hence lying on the past lightcone. Many objects, whose 
worldlines have not yet crossed our past lightcone, are not visible yet, so a large fraction of the Universe 
is unreachable to observation. 



Let us now turn to what such a model consists of. 



The now- standard model of cosmology is called after Alexander Friedmann, Georges 
Lemaitre, Howard Percy Robertson and Arthur Geoffrey Walker, who first introduced 
and discussed the corresponding metric and applied it to the Universe. The correspond- 
ing metric and space are called accordingly the FLRW metric, although for some unclear 
reason the "L" is often omitted... 



The FLRW model is mostly based on 4 basic assumptions: 

Theoretical prejudice and framework. Gravity is the leading force driving the dy- 
namics, and we shall describe it by means of the General Relativity (GR) theory of 
gravitation. Based on the equivalence principle, it is an extremely well tested theory 1 , in 



GR is accurate at the level of 10 , i.e. it compares with QED as far as predictions are concerned. This 
constraint is obtained by measuring the orbital period variation of the binary pulsar system which implies 



The FLRW model 



Theoretical hypothesis 



particular in the Solar System in which it serves as a reference for any alternative theory 
(scalar-tensor or MOND, for instance). 

Gravity it is the only known unscreened long-range force, and thus appears to be very 
well suited to describe the largest scales and even the Universe as a whole. In assuming 
GR to hold on these scales, we suppose the locally derived laws of physics apply and 
can be extrapolated. On the other hand, if anything were to go wrong in our description, 
that would lead to a natural testing ground for GR. 

The other interactions are assumed to be well described by the standard model of 
particle physics minimally coupled to gravity. This is achieved through the metric 
factor present in particular in the derivative terms: for instance, for a scalar field, the 
microscopic Lagrangian will contain a term of the form Xkinetic = -jg^d^dyc/). Hence, 
the fundamental action we shall be interested in reads 



It should be noticed again that any departure from this theoretical framework translates 
into observations different from the expectation, hence providing a way to test the 
validity of Eq. (4). In particular, scalar-tensor theories that would be equivalent to GR 
on Solar-System scales or for large cosmological times could originate (either in scale 
or time) very far from GR, and that could lead to observable consequences. 

Reasons for doubting the validity of GR in astrophysics and cosmology include the 
flat rotation curves of galaxies and the currently observed acceleration of the Universe. 
At least at a phenomenological level however, they can be described by GR providing 
extra "stuff" (dark matter and energy) is added to the matter content to which I now turn. 

Matter content. Once the theoretical framework is fixed, one needs to set the matter 
content, i.e. the right hand side of Einstein equation (1). Observations, made only over 
luminous matter (and hence not precluding a priori any dark component) on the past 
light cone, reveal a single class of objects, the luminous ones! Therefore, we need to 
model not only those observed objects, but also any other component that we would not 
actually be able to see. 

The typical distance scales involved are the galaxy characteristic size, of the order of 
10 6 light -years, and that of galaxy cluster, namely 10 8 light-years. Hence, we do expect 
some amount of clumsiness on scales of these orders: the large scale structure of the 
Universe, being supposedly insensitive to the small scale effects, will then be defined on 
scales larger than 10 8 light-years. On these scales, we will suppose the matter content 
to form a perfect fluid with normalized 4-velocity u 11 (g^ v u^u y = -1 with a metric with 
signature -2) and stress energy tensor 



where p(x, t) and p(x, t) are the energy density and pressure. The dynamics is usually 
imposed by setting V /U r /iV = 0, but this relates the time evolution of p and p in a contrived 



S = fd 4 X V=£ { TT^- (R - 2A) + Xmatter [<P(x), ^W, ■ ■ ■ , g^v] 




(4) 



T^ = (p + p)u^u v + pg tiV , 



(5) 



emission of gravitational waves in exact agreement with GR. This measurements led to its discoverers, 
Russell A. Hulse and Joseph H. Taylor, sharing the Nobel prize for physics in 1993. 



FIGURE 2. We observe an isotropic distribution of matter, and that can correspond either to an homoge- 
neous distribution (left) or to a spherical shell-like structure centered on our location (right). Copernican 
principle demands the first option, but such a philosophical posture, however well justified, needs be 
verified experimentally; this can nowadays be done by 3 dimensional observations using redshift data. 



way: one needs to impose another relation, called the equation of state, expressing the 
pressure as a unique function of the energy density. In practice, a few simple cases are 
considered, always assuming a linear relationship, i.e. p = wp, with w a constant called 
the equation of state parameter. 

Observations reveal the typical relative velocity between galaxies (the point particles 
in the fluid element description) to be of the order (v ga i) - 200km • s" 1 ~ 10 -3 in units of 
the speed of light. Therefore, the mean kinetic energy relative to the mass can be eval- 
uated as ^kin/p - ^( v g a i) ~ 10~ 6 , and this also provides a measure of the numerically 

expected value of the ratio between pressure and density: p/p ~ ^(v?^) ~ 10 -6 . There- 
fore, the fluid made up with the galaxies and any similar behaving fluid (dark matter) 
will be described by a pressureless gas, i.e. w m ~ 0. 

We also observed that the Universe is filled with some amount of radiation, whose 
stress energy tensor is traceless, thus implying w T = ^. Finally, a cosmological constant 
term can be described by writing = -(A/8^G N )^ y , and a direct comparison with 
(5) then shows that this implies pa = ~Pa, m other words wa = -1. Amazingly, this 
extremely simple set of 3 constant equation of state fluids suffices to describe the 
evolution of the Universe for the previous 13.7 billions of years with percent accuracy! 

Symmetries. Without symmetry assumptions, it is impossible to solve the full GR 
equations, even with a given (and simple) stress energy tensor such as that presented 
above, and so one needs to make even more simplifying assumptions, again based on 
observations. Those reveal the distribution of matter and radiation to be essentially the 
same in all directions. In other words, we see a space which appears isotropic. Figure 
2 then implies at least two options following from these observations, of which the 
simplest is homogeneity (but spherical symmetry has also been studied), to which I 
will stick for now on. It should be emphasized that both homogeneity and isotropy are 
concepts whose validity in cosmology makes only statistical sense, and it is in this sense 
that they must be verified whenever possible. 



At the level of Newton classical theory, homogeneity, stating that each point of 
space is similar to any other at each instant of time, is well defined. In GR however, 
the previous sentence is absolutely meaningless, and requires that a 3 + 1 (space and 
time) slicing is done, hence generating a one parameter (time t) family of spacelike 
hypersurfaces Ej. Homogeneity then is rephrased by saying that for any two points in 
S f , there exists an isometry taking one to the other. Isotropy on the other hand states that 
at each spacetime event, an observer moving with the cosmic fluid (comoving observer) 
cannot distinguish one direction of space from another one. One sees that the two notions 
are quite intricate, even though one describes a property of spatial hypersurfaces, while 
the other involves time development; this is due to the nature of our observations, always 
done along a light cone, hence mixing space and time measurements. Figure 3 clarifies 
these statements. 
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FIGURE 3. Homogeneity (left) and isotropy (right) can be given exact mathematical meaning in GR: 
any two points P and Q can be related in an invariant way through an isometry in the hypersurfaces 
Ej (homogeneity — > generalization of translations in space), and for any point P and two spacelike 
orthonormal vectors e x and such that e^u^ = 0, there exists an isometry transforming e\ <-> ei (isotropy 
— > generalization of rotations). 



Assuming homogeneity and isotropy means that h(t) = g\ Pei: , restriction of the full 
metric g to the hypersurface 2 f , only depends on time t, so that l, t is a 3 dimensional 
homogeneous and isotropic space with induced metric 
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as can be checked by direct calculation of a vector lying in E f or orthogonal to it. 

Let ^R^vap be the Riemann tensor associated with the metric h^ v on S r . By its 
definition (2), it is symmetric in the exchange of the pairs of indices {ftv} and {a/3} and 
so can be seen as a map L of the vector space of 2-forms onto itself: setting A = {/uv} and 
B = {a/3}, the matrix L£ representing the 3 dimensional Riemann tensor is a symmetric 
matrix and thus diagonalizable. If its eigenvalues were not all equal, then one of them 
could be used to generate a privileged direction, in contradiction with the hypothesis 
of isotropy. Hence, we have L = Kl, where K 6 R can only depend on time, and 1 the 
identity in the relevant space. 

Moving back to 4 dimensional indices, we can write the 3 dimensional Riemann 
tensor in the form 

(3) Rfiva/3 = Kit) (h^hyp - hayhrf) . (6) 



Let us see the meaning of this expansion for K > to begin with, and embed the 3 
dimensional space in a 4 dimensional Euclidian space with coordinates x, y z and w. 
A constant positive curvature space is a 3-sphere of radius a whose point locations are 
given by 

x z +y 2 +z Z + w 2 = a 2 . (7) 
In spherical coordinates defined by 



x = acosx, 

y = asin^cos^, 

z = asin^sin#cos</?, 

w = asin^sin^sin^, 



(8) 



differentiation of Eq. (7) then provides the 3 dimensional metric in the form 

d (3) 5 2 = (dx 2 + dy 2 + dz 2 + dw 2 ) | = a 2 (t) [d X 2 + sin 2 ^ (d<9 2 + sin 2 6dip 2 )] , (9) 

where in the last line we have put back the possible time dependence of the overall 
spatially constant curvature. 

Similar considerations with negative (3-hyperboloid) or flat (Euclidian) space permit 
to rewrite the overall 4 dimensional metric in the special FLRW (at last!) form 

ds 2 = g^vdx^ dx y = (h^ v - u^iiv) d^d^: v = - (u^dx^ + h lIv dx fI dx y = -dt 2 + a 2 (t)y t jdx l dx j , 

(10) 

where the spacelike part of the metric is 

7ij dx l dxj = d X 2 + fi(x)d0 2 , with /^=«7C- 1/2 sin(V^), 

where d£l 2 = d6 2 + sin 2 8d(f 2 is the usual solid angle element and the function fix is to be 
continued for vanishing (lirn^o 7ft = X) or negative [fix<o —* (-7<")~ 1/2 sinh( V-^C^j] 
values of < K. 

In the above relations, we have written 'K to distinguish from the function K(t) giving 
the 3 dimensional Riemann tensor. In fact, it is always possible to renormalize the spatial 
coordinates in such a way that the scale factor a(t) has the dimension of length, so that 
'K can take one of the possible values 'K e {0, + 1}. This is the choice we will assume for 
now on. 

To finish this paragraph, I suggest to the reader to try and show, as an exercise, that 
the spatial metric can be cast in the equivalent forms 

. i, i dr2 2,^2 dP + PdQ 2 

7ijdxdx = + rda = 1 P {U) 

l-Kr (l+7C4^ 2 ) 
by means of changes of coordinates x —* r —> £ to be determined. 



Topology. GR is a local theory which thus says nothing about the global structure of 
the Universe. String theory, of which GR is a low energy limit, teaches us that some 
dimensions may be compact, and in fact need be so in order for space to appear 3 
dimensional on the scales available to experiments. Therefore, in principle at least, 
it is possible that the large dimensions we happen to live in could also be compact, 
leading to a non trivial topology. Although this has been studied in details, I shall not 
consider any further this hypothesis for at least two reasons, one experimental and the 
other theoretical. 

First, there is no data, to date, that would induce us to think a large scale non trivial 
topology is needed. Of course, some compact models can be made compatible with 
the data, or even improve the fit, but they are degenerate with other models and the 
improvement is not really statistically significant. 

Second, an argument in favor of compact large dimensions could be to invoke com- 
pactness for all dimensions; in this case, the expected phase of inflation would make the 
large dimensions much larger than the current Hubble scale (size of the observable uni- 
verse), unless a disturbingly severe fine-tuning is applied. If a non inflationary scenario 
is implemented, then another fine-tuning is necessary in order to explain why the lattice 
size of the compact dimension should be, today, of the order of the Hubble scale (only 
case not yet ruled out by the data but still leading to observable predictions). 

Having settled the framework, let me move on to the dynamics of our Universe. 



The dynamical Universe 

The framework developed above permits to write down explicitly the Einstein equa- 
tions as a set of relations between a very small subset of dynamical quantities, namely 
the scale factor a(t) and the density of the fluid p(t). In order to derive these equations, 
we need to calculate all the relevant geometrical quantities. 



Geometrical quantities 

The Einstein equations involve in a non trivial way the Riemann tensor and its 
byproducts, namely the Ricci tensor and scalar and the Einstein tensor itself. Those 
are all built from the metric connections and ultimately from the metric itself. It turns 
out that the cosmic time t introduced earlier is not the most convenient time parameter, 
especially when the spatial sections are flat (which is observationally the case), and we 
usually introduce a dimensionless time, called the conformal time rj as it renders the 
metric conformally flat. It is related to the cosmic time by 

adrj = dt => ds 2 = a 2 {if) {-drj 2 + y t jdx'dx J ) = a 2 {if) (-drj 2 + dx 2 + dy 2 + dz 2 ) , 

v ' K-tO v ' 

(12) 

where in the last stage we have taken the limit *K —* to make the Minkowski metric 
apparent. 



To simplify matters, we define derivatives with respect to times as / = df/dt = 
df Hadif) = a~ l f. Then, setting H = a/a and *H = a' /a = a = aH, we obtain the only 
non vanishing connection coefficients as 

Y\ r a 2 Hy lh V tj = H6) and Y) k = y) k (13) 

in cosmic time, and 

T%,=<H, Tl = <H 7ij and ^ = ^6), (14) 
in conformal time. From these, one derives the non vanishing Ricci tensor components 

R„ = -3-, R ti = and = a 2 7 J- + 2H 2 + 2% ), (15) 
a J \a a 1 ) 

leading to the Ricci scalar R = 6^H 2 + a/a + ^K/a 2 ^. 

Combining these, we finally obtain the Einstein tensor as 

°:=- 3 (« 2+ f). °:=° a » d G i=-4(^+» 2 +|)- < 16 > 

This provides the left hand side of Einstein equations (1). 



Friedmann equations and the cosmological parameters 

With the geometric quantities derived for the FLRW metric, and the stress energy 
tensor (5) restated in matrix form as T^ y = diag(-p, p,p,p), it now remains to equal it 
to (16) to obtain the Friedmann equation, which the reader will straightforwardly check 
they can be cast in the form 

tf = !^,-*4 (17) 

3 a 1 3 



for the constraint, and 



4ttG n . . . A 



, (P + 3 P ) + - (18) 
a 3 3 

for the dynamics. Deriving Eq. (17) with respect to time, taking into account Eq. (18) 
and reshuffling the various terms involved yields the fluid conservation equation 

p + 3H(p + p) = <= V^ v = 0, (19) 

as expected since the latter conservation is not independent of the Einstein equations 
from which (17) and (18) stem. 

In terms of conformal time, the previous set of equations read 



p' + 3'H(p + p) = 0, 



(20) 



for the conservation equation, 



-K 2 + 7C= — -^p + -\a 2 (21) 



for the constraint, and finally 



4^G N A 
■^(P + 3 P ) + 3 



a 2 . (22) 



There exists a special solution, which happens to be realized in our Universe, at least 
so seem to say the data, namely that for which the spatial curvature 'K vanishes. It defines 
a density, called the critical density p c given by 

8?rG N p c 

in terms of which one can express all densities in a dimensionless way. For each fluid 
component but the cosmological constant, one can set Q. a = &7rG N p a /(3H 2 ) = p fl /p c ; 
we also introduce an equivalent curvature "density" as Q«k = - < K/(a 2 H 2 ) and finally 
£2a = A/(3// 2 ), and then the Friedmann constraint simply reads: 

^fi fl + Q A + % = l, (24) 

a 

so the Friedmann equation is understandable as an energy budget: all possible contri- 
butions basically sum up to 100%! Numerically, the Hubble constant today is measured 
to be of the order of H = 100/zkm • s" 1 • Mpc" 1 , where h = 0.704 + 0.025. Similarly, 
the relative densities are also measured in units of the critical density, estimated as 
p c - 1.9 x 10 _29 g • cm -3 ; they frequently are found expressed as p® = C$h 2 to account 
for the indeterminacy of the Hubble expansion rate as well as on the density parameter 
itself, the subscript "0" meaning the present-day value. 



Special solution: matter and radiation 

With a varying equation of state w(i) and a scale factor a(t), which is a monotonic 
function of time, it is always possible to parameterize all functions of time as functions 
of a, and in particular w. On can then formally integrate the conservation equation as 



p[fl(0]=Piniexp|-3 J [l+w(a)]dlna^ = p M ( — | , (25) 



■3(1 +w) 

w^cst' ~~\a; r 



which gives an exact solution for the constant equation of state situation. This is pre- 
cisely the case when matter (w — > w m = 0) or radiation (w — » w r = ^) dominates over 
everything else. Eq. (25) then shows that matter scales as p m oc a~ 3 , as expected from 
mass conservation in an expanding volume, while radiation gets an extra power, scaling 



as p r oc a" 4 , due to the redshift of its wavelength. Now consider an initial condition con- 
sisting of given relative amounts of matter and radiation. When the Universe begins its 
evolution, with a small value of the scale factor, radiation dominates and the total den- 
sity is p tot ~ p r until it gets caught up by the dustlike matter. This remarkably accurate 
picture for the Universe density evolution is illustrated in figure 4. 
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FIGURE 4. Top - Evolution of densities: the Universe begins dominated by radiation, whose density 
decreases faster than that of matter, so the latter ultimately dominates. Not shown is the final phase of 
domination by a cosmological constant which, as its name indicates, behaves as a constant. The point 
at which radiation and matter contribute equally is, not surprisingly, called equality. Bottom - On the 
same scale, matter density is depicted together with a typical light ray, whose mean free path is initially 
much shorter than the Hubble scale, as e.g. during nucleosynthesis; as the matter density gets smaller 
and smaller, the mean free path eventually becomes larger than the Hubble scale after what is therefore 
denoted decoupling. The Universe becomes transparent to this radiation we now observe as the microwave 
background. 

The meaning of the equation of state is clarified when one considers a perturbation 
propagating in the fluid. As is well known in fluid dynamics and as we shall also discuss 
later, the sound velocity c s is given 2 by cl = dp /dp = p' /p'. It can be shown (and the 
reader is encouraged to do so!), that the relation 

W = -37/(1 +w)(cj-w) ( 26 ) 

holds, so that a constant equation of state means w = cl . 

With the solution for the density as a function of the scale factor and the equation of 
state given, it is an easy matter to solve the Friedman equation. For a vanishing spatial 



In fact, it should be partial derivative for constant entropy. 



curvature "K = 0, one finds that if w ^ - 1 , the solution goes like 

flOC? 2/[3(l + w) ]oc?/ 2/(l + 3w) ^ a.oc^tocrj and a m oc t 2/3 oc f , (27) 

where I emphasized the particular pressureless dust and radiation dominated solutions. 
In the special case of a cosmological constant with w = -1, one finds 

p = =^ H = cst =^ aoce Ht oc—, (28) 

and one has an exponentially accelerated expansion; note in that case, which will later 
correspond to the inflationary situation, that the conformal time is negative, with the end 
of inflation being for the limit when rj — » 0. 

Limitations of the standard model 

The model developed above gives a quite accurate description of the history of the 
Universe, but its success actually raises a few questions that find no answer in its own 
framework. 



Puzzles 

Singularity. The first troubling issue is also the only one that has, in the inflationary 
paradigm, not received any answer, namely the fact that whatever solution of Einstein 
equations one comes up with that fits the available observational data does begin with 
a primordial singularity: at some point in the past, there is always a time t s i ng at which 
<2feng) — » 0, meaning all the geometrical tensors diverge, so the theory itself simply does 
not make sense anymore! One can however argue that GR is not designed to handle 
extremely high energies so that a cutoff, at the string or Planck scale, should be applied, 
above which the theory will (wishful thinking) be regular. 

Horizon. The question of the horizon is more involved in a way, as no hand-waving 
argument can be similarly invoked to cure it. It relies on the observed fact that light 
emitted at decoupling (see figure 4) is homogeneous up to 10~ 5 . Although this looks 
like a mere consequence of the cosmological principle, it is actually weird because of 
the previously discussed singularity problem: the existence of a primordial singularity 
implies a Big-Bang, i.e. a point in time at which the Universe expansion starts, so that 
there was a finite amount of time for a priori initially causally disconnected regions to 
thermalize. When one estimates the number of such regions, one finds some 10 5 of those 
at decoupling, implying a predicted isotropy over angular scales smaller than roughly 
one degree on the sky only! Figure 5 illustrates the issue. 

Flatness. Finally, the flatness problem is based on the fact that the observed flat 
spatial section (2 a ^« + = 1, i-e- Oft = 0) is actually an unstable fixed point: in the 
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FIGURE 5. The horizon problem. Left - spacetime diagram beginning at the Big-Bang (i.e. the sin- 
gularity). Light emitted in A and B in all possible directions reach regions at the last scattering surface 
(surface of decoupling when light stopped scattering and started propagating unaffected) that never were 
in causal contact. Yet, they appear to have exactly the same physical properties. Right - Angular repre- 
sentation of the same thing: the Big-Bang singularity is now represented from our point of view by the 
infinite redshift sphere. Calculating the horizon size at decoupling gives one degree on the sky, which is 
thus the maximal angular scale over which one might expect to measure an isotropic distribution. 



absence of a dominating cosmological constant, deriving Eq. (24) with respect to the 
scale factor yields 

^ = (3w+l)(l-0 K )a K , (29) 
dlna 

whose solution, for a constant equation of state, reads 



Q<K (^obs) = &x 



^obs 



l+3w 



-1 



(30) 



where £2™ = Q.<x(a m i)- In order to observe now Q^-(a Q bs) ~ 0.1, one then needs to 

demand that at equality (a bs/«eq - 10 4 ), |OS| £ 3 x 10 -5 . which represents already a 
quite substantial amount of fine tuning if Q.<x is to be an arbitrary initial condition. It 
becomes even worse of course if one assumes initial conditions ought to be imposed at 
one Planck time after the Big-Bang singularity itself, as the requirement then becomes 
l^pianckj < iQ-60. fa[ s j s an unacceptably large amount of fine tuning! 



Categories of solutions 

There are nowadays two categories of solutions as far as I can tell, one widely 
accepted and usually set as part of the standard cosmological paradigm, namely inflation, 
and a contender based on a contracting phase and a bounce. Inflation does not address the 
singularity question, while a bounce is in danger of producing too much shear during the 



contraction. Moreover, inflation can be easily implemented using a simple scalar field, 
a de Sitter phase actually being an attractor in the equations of motion, while a bounce 
is almost impossible to implement in the framework of 4 dimensional GR unless the 
spatial curvature is positive. This makes inflation more appealing to most people. 

Inflation. The Flatness problem originates from Eq. (29) and the fact that Q,<x = 
is an unstable fixed point for this equation. In fact, this is only true provided w > 
If this condition is not fulfilled, as for instance is the case of a cosmological 
constant domination, then £1% = becomes instead an attractor. So it suffices to include 
a sufficiently long phase during which Q.<x —* 0, then followed by the usual radiation and 
matter domination, to keep Q<x close to zero even after a long time of regular expansion. 
What is the meaning of this solution? 

Equation (18) in the absence of a cosmological constant shows that if w < —h, i.e. 

if p < —\p, then a changes sign and the expansion is accelerated. This is why this 
solution was called inflation. Very often, it is implemented by means of a "slowly 
rolling" scalar field, dubbed inSaton, i.e. a scalar field whose dynamics is dominated 
by the potential term, naturally leading to w -1. As a result, inflation is achieved by 
an almost exponential growth of the scale factor. 

Having a phase of accelerated expansion actually also solves without any further 
assumption the horizon problem. Indeed, the horizon size is a global quantity whose 
definition involves the overall history of the Universe through 



where t m i is the origin of times. Note that for a power-law expansion such as during 
most of the history of the Universe [see Eq. (27)], i.e. if a oc t a , the horizon scales as 
dft = t/(a+ 1), which is then roughly the same as the Hubble expansion rate H~ l = 
a/a = t/a. This, plus the fact that "Hubble" and "Horizon" begin with the same letter, 
has led to a confusion in many works between the two quantities. I will come back to 
that point later. 

When an almost exponential phase of inflation takes place, the Hubble radius is 
roughly constant, while the scale factor grows exponentially. The horizon size is then 



where the last inequality assumes t » ti n \. What happens then is that the horizon size 
grows much faster than the Hubble scale so that all scales end up having time to be in 
causal contact. 

To be complete with the inflationary scenario, apart from its prediction of an almost 
scale-invariant spectrum of primordial fluctuations, I should like to mention that it is also 
the only known way of naturally reducing any initial amount of anisotropy. However, 
in order to reach an FLRW Universe, one also needs to impose a sufficiently smooth 
initial patch, i.e. even though inflation substantially alleviates the question of initial 
inhomogeneities, it does not actually answer it. 




(3D 



rf H = ife H(f - f - ) -ll» J f/- 1 , 



(32) 



Bouncing scenarios. Eq. (29) can also be rewritten as 



-77- = -273, (33) 
at a i 

emphasizing once again that a phase of accelerated (a. > 0) expansion (a > 0) will drive 
£1<X to vanishingly small values. Another way is possible, consisting in reverting all the 
signs of the previous argument, therefore using a decelerated (a < 0) phase of contraction 
(a < 0) ! Since we observe the current phase to be expanding, this implies that a transition 
between H < and H > took place, a bounce. 

One might however argue that, as I said earlier, a positive spatial curvature is required 
to implement such a bounce in GR. One quick answer to this argument is that GR 
may not be valid at that time... after all, inflation also requires some extension of GR 
to account for the primordial singularity. Note in passing that the singularity is easily 
gotten rid of in the bouncing scenario since time can be pushed back as far as one wants, 
in principle even to infinity. This actually also solve the horizon problem, since it can 
easily be made infinite. To anyone not willing to extend GR in any way, one could also 
argue that it suffices to have a very long contraction phase during which Q.<x —> 0, and 
manage that the bounce itself is not very asymmetric, so that even though Q.<x can grow 
large during the bounce, it will recover a value after the bounce that is not very different 
from the one it had before. 

The bouncing scenario however is plagued with an anisotropy problem: if one con- 
siders a initial shear, however tiny, it will grow very large during either the contracting 
phase or during the bounce itself. This is at least true in simple models, but more sophis- 
ticated scenarios have been proposed that tame these unwanted large growth. 



Structure formation: perturbation theory. To obtain a complete description of the 
Universe, one would, at this point, need to include thermodynamical evolution of all 
the relevant quantities, taking into account interactions to describe in a reasonable way 
the phases of nucleosynthesis for instance. I shall not embark in this direction, and will 
instead concentrate on the perturbations over this background: those in fact provide a 
bonus for the inflationary scenario, as by demanding the inflaton to be in vacuum and 
allowing it to have quantum fluctuations, the ensuing evolution transforms the Universe 
in a particle producer, those particles then later behaving as large scale fluctuations 
seeding the formation of structures. 



PERTURBATION THEORY: SVT AND THE GAUGE ISSUE 

It is widely believed that large scale structures formed out of primordial seeds upon 
which gravitational collapse acted to produce dense objects. This theory is quite well 
verified, as numerical simulations starting with an initial over-density function satisfy- 
ing scale-invariant statistical properties manage to reproduce the statistical properties 
observed in large scale structure surveys. It lacks however a crucial ingredient: what is 
the seed origin? 



Introductory remarks: the Jeans length 
and Newtonian perturbation theory 



Newtonian physics allows to understand the origin of gravitational collapse in the 
expanding Universe in a phenomenological way: one simply assumes that Newtonian 
gravity holds, but also that the Universe is expanding, i.e. that the actual distance r 
between objects increases with time. One then has r = a(t)x, where x is the relative 
position of the object in a local coordinate system and a(t) the scale factor discussed in 
the previous section. The total velocity then consists in two pieces, 

— = v = ax + a— = aHx + u , (34) 
ot ot 

background peculiar 

where the first term represents the background cosmic flow involving the Hubble rate 
H, and the second the peculiar velocity, i.e. a relative velocity that one can treat as a 
perturbation. 

Similarly, the density field is expanded as 

p(x,t)=p(t)[l+6(x,m, (35) 

and the continuity equation becomes 

+V r -(pv) = => + 3Hp+^V x (pu) = 0, (36) 

where the '3//' term comes from changing the coordinate r to x. To zeroth order, Eq. 
(36) implies p + 3Hp = 0, which merely reflects that matter scales as p oc a" 3 , while the 
first order yields 

(5+-V-[(1+<5)h] = 0, (37) 
a 

where for now on we assume all spacelike derivatives are with respect to the 'comoving' 
coordinates x; I shall accordingly subsequently omit the index x. 
Combining Eq. (37) with 

^ + //« + - (H • V)m + - 1 -VP + V<D| = 0. (38) 
ot a a\p ) 

which is nothing but the Euler equation for a fluid with pressure P in a gravitational 
potential O (satisfying the Poisson equation AO = AnG N p8) in comoving coordinates, 
and linearizing, one finds 

6+ 2H6^ - -|A5 =4nG N p6, (39) 

expansion ^avitT^ 
pressure b J 



showing the Newtonian evolution involves three distinct effects, namely the damping 
of any perturbation due to the expansion, the propagation of sound waves due to the 



pressure terms, and finally gravity itself. In Eq. (39), the sound velocity c s is defined as 
before through c 2 = (dp I dp) s where the fluid entropy S is held constant. 

Expanding in Fourier modes (A —* -k 2 ) and defining the physical wavenumber k p = 
k/a, one obtains, forgetting for the moment the expansion (i.e. setting a —* C te ) 



6 + (c 2 s k 2 p -4nG N p)6 = 



tfocexp 



4nG N p 1 



A 2 



(40) 



where wavelengths are defined by A = 2n/k, and A j = c s sJn/(G N p) is the celebrated 
Jean's length separating regimes of wave propagation and gravitational instability: for 
long wavelength, A > A v the density is growing exponentially with time, signaling a 
collapse, while for small wavelengths A < A v the density oscillates as the sound wave 
propagates smoothly. Taking into account the overall expansion does not change this 
result qualitatively, it merely changes the functional dependence of the density with 
time, not the fact that there is a regime of unlimited growth and another of oscillations. 

Having settled the stage, let me now move to the real issue, namely that of GR 
perturbations in FLRW Universe. 



3+1 decomposition 

For now on, I will mostly consider the conformal time 77, in terms of which the 
subsequent exposition is probably clearer. Then, the Friedmann (Einstein) equations are 
given by (21) and (22). Having completely fixed the background, we can now move on 
and expand around this background. 



Perturbative expansion 

Our starting point is the action (4) or, in practice, Einstein equations (1). As we did 
obtain the homogeneous and isotropic solution, we write it as gjjirf), leading to the 
corresponding Einstein tensor G^J (rf), itself sourced by the background stress-energy 
tensor T™ {rf). We then write the full metric as 

4v U 07. *) = C 07) + eg™ (17, X) + l -e 2 g^ (17, *) + ■■■ , (41) 

where the dots contain all higher order terms. We then assume that e is a small parameter, 
as data indicate it to be the case on sufficiently large scales, i.e. on scales larger than 
roughly 200 Mpc. With Eq. (41) and the definition of the Einstein tensor, one can express 
it in the same way, namely 

Gjf (77, x) = G$ (77) + sG$ (77, x) + l -s 2 G^ (77, x) + • • • , (42) 



where G^ y (77) is given by (16). 

Similarly, we expand the stress-energy tensor (5) as 

T™ 1 (77, x) = T$ (77) + eT™ (77, x) + X -s 2 T§ (17, *) + ■■■ , (43) 

which amounts to expanding p, p and the fluid vector u^. Providing the series in powers 
of s makes sense, it now suffices to expand both sides of Einstein equations and identify 
the terms, order by order. 

In practice, there is no s parameter, and we merely expand all relevant quantities as 
"background" + "something small" which we then calculate. The metric g^ v = g^v + dg^v 
(denoting for now on the background quantities by an overbar) will read 

ds 2 = a 2 (77) [- ( 1 + 2A) drj 2 + IB^dx 1 + (y y + fcy) dx i dx s ] , (44) 

whose "e —> 0" limit would give (12) back. Note that the quantity g^ y = g^ y + dg^ 
should be the inverse of the above metric, so that demanding g^gva = 5^ a , we obtain 

The stress-energy tensor (5) has background values obtained with the choice = 
a~ 1 6 t ^, i.e. = —ctffy, and we also demand that the timelike vector = uF + 6u^ 

be normalized at all orders, leading to 6u? = a~ l (-A,v^, thus defining v', and du^ = 
a(-A,Vi + Bi); we see that it depends on the metric perturbation. 
Gathering all terms for the stress-energy tensor, we finally obtain 

8T m = a 2 p (6 + 2A) , 6T ¥ = -a 2 p [(l+w) v,- + B{\ and STy = a 2 p | ^y (7 + hijj , (45) 

where in the last term we have omitted a possible anisotropic stress contribution. The 
equation of state itself is perturbed assuming now that the pressure is a thermodynamical 
function of both the energy density and the entropy, if any: p = p(p,S). The pressure 
perturbation then reads 

5p = cfa + t6S = c 2 5p + pY = c 2 & 5p + 6p nad , (46) 

where I indicate the most frequently used notations. The last one refers explicitly to 
the "non adiabatic" component of the pressure, which is proportional to the entropy 
variation SS . 



Scalar, vectors and tensor components 

In the perturbative expansion, we see appearing ordinary functions and indexed ob- 
jects. The former transforms as scalars on the spatial hypersurfaces (recall we have an 
explicit 3+1 decomposition), while the latter transform as either vectors or rank-2 ten- 
sors of the spatial sections. 

Making use of the covariant derivative associated with the metric y ;j , which we call 
Di, one can decompose all relevant quantities in terms of pure scalar, vector and tensor 



modes: for instance, the vector 5 ; appearing in the metric (44) can always be written as 

B i = D i B + B i , where D'fl, = 0, (47) 

thus exhibiting a scalar function B and two divergenceless vector degrees of freedom 5,, 
recovering the three initial vector degrees of freedom. In a more common hydrodynami- 
cal framework for instance, that would be equivalent to splitting the velocity field v; into 
a velocity potential D ; v and a vorticity term v,-. 

The same technique applies to the tensor quantity hij which we write as 3 

h ij = 2[Cy i j + D (i D j) E + D (i E f) + E ij \ with D i E j = and D i E ij = = E i j , (48) 

where now the tensor Ejj is not only divergenceless but also traceless. This way, the 
10 degrees of freedom of the metric are now split into four scalars (A, B, C and E), 
2 vectors (5 ; and £,) of 2 degrees of freedom each, and one tensor Eij, also having 2 
independent degrees of freedom. The main interest of this Scalar- Vector- Tensor (SVT) 
decomposition is that, at linear order, they all decouple, and one can thus treat the scalar, 
vector and tensor modes independently. 



The gauge issue 

GR is diffeomorphism invariant, i.e. it is constructed in such a way that general 
coordinate transformations leave the equations unchanged. This implies that out of 
the 10 degrees of freedom discussed above, 4 are essentially irrelevant as they can be 
absorbed into a coordinate transformation. When applied to the special background + 
perturbations case, this invariance is no longer an actual coordinate transformations since 
the background is kept fixed; it is then called a gauge transformation. Let us see in more 
details how it works. 



Metric fluctuations 

Suppose I change the coordinates x^ to a set of new coordinates jc^ related with the 
previous ones by an infinitesimal translation, i.e. i-» x^ = xF + ^, where ^ are small 
quantities. General covariance then implies that the equations of motion have the same 
form when expressed in the "new" coordinates x? or the "old" ones x^. In particular, the 
line element, namely ds 2 , should have the same structure under a gauge transformation. 
Therefore, we set 

ds 2 = a 2 (r]) [- ( 1 + 2A) dfj 2 + IBidfjdx 1 + (y (7 + h t j) dx ! 'djc y ] , (49) 

and by gauge invariance, we require that ds 2 = ds 2 , after having SVT-decomposed 
the transformation through f] = rj + T and x l = x' + D'L + U. We find the following 



3 We denote by round parenthesis the symmetrized part of the relevant tensor, i.e. fyj) = \ [fij + fjj 



transformation laws: 



A=A-(T' +<HT), B = B-(L'-T), C = C-<HT and E = E-L (50) 
for the scalar quantities, 

fy = B i_li> and & = (51) 

for the vectors and finally E[j = Ejj. The last identity could have been obtained without 
any calculation from the vectorial nature of the transformation: tensor modes, also called 
gravitational waves, are naturally gauge invariant. 

Equation (51) can easily be reshuffled into E 1 ' - B l = E 1 ' - B 1 , so that the quantity 
O' = E" - B l is gauge invariant. On the scalar side, similarly, one finds that (^B-E'^j = 

(B - E') + T, so that [A + (fl - E')] = [A + (B- E')] -<HT, and finally that the quantity 

d>=A + (B-E')' +<H{B-E') (52) 

is also gauge invariant. I leave it as an exercise to show that 

y ¥ = -C-'H(B-E') (53) 

closes the set of gauge-invariant variables consisting of two scalars O and called the 
Bardeen potentials, two vectors O 1 and two tensors for a total of 6 gauge-invariant 
quantities, as expected from the original 10 quantities and 4 possible gauge choices. 



Choosing a gauge 

One can play the same game with the stress-energy tensor and obtain transformation 
rules by expressing it in one frame or the other through the usual transformation rule of 
a rank-2 tensor. One finds 

Sp = Sp+p'T, v = v-L', v ! ' = v'-Z/' and 6p = 6p + p'T, (54) 

leading here also to a set of gauge-invariant variables 

6p N =6p+p'(B-E'), 6p N =6p + p'(B-E'), V = v + E' and V*' = v' + B\ (55) 

given here an only one example of such a combination. 

From that point on, one can write down Einstein equations and solve them: just like 
in electromagnetism, one merely needs to fix a gauge. There are many gauges that have 
been used in the literature, and I list a few of them here. The first I want to list shows that 
the gauge-fixing choice is, just like in electromagnetism again, not necessarily enough: 
it is the so-called synchronous gauge, in which only spatial sections are perturbed. In 
other words, it is defined by assuming that the proper time of a comoving observer is 
cosmic time, and this translates into setting A = and B[ = 0. Because of its definition, it 



is a quite intuitive gauge, but it is not completely fixed: setting r) = f(rj) or x l = /' (xA for 

arbitrary fonctions / and /', one remains in this gauge (A = and 5 ; = are still valid). 
This leads to possibly spurious solutions, and hence to mistakes! 

Another frequently used gauge in the case of a single fluid is one which follows the 
fluid's motion, so that one demands 6T. = 0. This is an interesting choice which becomes 
unfortunately ambiguous as soon as more than one fluid is involved. In this gauge, the 
variables 

6p c = 6p + p'(v + B) and 6p c = 6p + p (v + B) , (56) 

are the natural fluid variables to use. 

Another physically interesting choice is that which consists in demanding the curva- 
ture perturbation of spatial section to vanish, which amounts to setting C = E = and 
Ei = 0, so the quantities 

S/ f=Sp-p'^ and Sp^Sp-p'^, (57) 

reduce to their original values: these gauge-invariant variables are thus the density and 
pressure perturbations in the flat-slicing gauge. 

Finally, it seems also appropriate to use directly a set of physically relevant variables 
like those defined above, namely the gauge-invariant ones. The simplest way to do that 
is to impose the so-called longitudinal, or Newtonian, gauge, i.e. that in which the scalar 
part of gpy is diagonal so that we set E = B = 0. In this gauge, the potential (53) is 
the Newtonian potential. There are two possibilities to get to this gauge: one can either 
set E = B = from the outset (easy way) or work out all the equations and express all 
of them only in terms of the gauge-invariant variables (52) and (53). They both give 
the exact same results, of course. One sees that in this gauge, the density and pressure 
perturbations reduce naturally to those defined in (55). 



Perturbed Einstein equations 

We now are in a position to write down explicitly the Einstein equations to first 
order of perturbations in a meaningful way. The equations in the Newtonian gauge only 
involve gauge-invariant quantities, and I shall therefore restrict attention to those in what 
follows. Since the following section is dedicated to tensor modes, I will simply forget 
about them until then (remember they decouple at linear order anyway). 

The next-to-simple case is that of vector modes. In most cosmologically relevant 
situations, there is no anisotropic stress (fti = 0), so that the equations of motion of the 
vector modes are not sourced by anything. They take the form 

(A + 27C) <!>,■ = ^pa 2 (1 + w) V u (58) 

and, more importantly 

o; + 2Wi = pa 1 fit -> 0, (59) 



leading to the exact solution 6, oc a~ 2 , and consequently, thanks to (58), that V,- oc a 3w ~ l . 
It is a well-known (observational) fact that vector modes were negligible at the time 
of nucleosynthesis, so we may confidently set ||0 ; || 1 at z nuc i ~ 3 x 10 8 . This implies 
||Oj|| «; 10~ 17 now: apart in very special situations such as a contracting universe in a 
bouncing scenario, one can set the vector perturbation to zero. I shall not consider them 
anymore in what follows. 

We are thus left with scalar modes. Since those have been driving the gravitational 
collapse leading to large- scale- structure formation, they are definitely the most relevant 
modes to study and, indeed, they play the first role in most of the literature on the subject. 
Their time development is obtained through two independent sets of equations, the first 
relating density to pressure perturbations, i.e. Eq. (46), the rest being given by Einstein 
equations, the spatial part of which, proportional to 8T l . oc 8 l . for a perfect fluid, yielding 

yWjDji^-W) = 0: under the reasonable assumption that the scalar perturbations do 
not diverge at spatial infinity, this relation implies that the only possibility is to have 
*F = O, a condition which I will take as valid for now on. 

For the scalar modes in the longitudinal gauge, Einstein equations then read 

AO-3W -3( < H 2 - < K)® = 4nG N a 2 Sp N , (60) 
A(0' +W) = -4nG N a 2 (p + p)ViV, (61) 
O" + 3W + {m' + <H 2 - 7C) O = 4nG N a 2 6p N . (62) 

Equation (60) can be reformulated as (A + 37<") O = AnG n a 2 8p c , using (56) and (61). This 
Poisson equation (up to the spatial curvature term) shows that the Bardeen potential is 
essentially the ordinary Newton potential if the density perturbation is expressed in the 
comoving gauge. As we shall see later, the matter perturbations in the different gauges on 
scales smaller than the Hubble radius are basically the same, so the sub-Hubble Bardeen 
potential indeed reduces to the Newtonian one (hence the notation O). 

Now, using Eq. (46) to express the pressure perturbation in terms of the density, and 
then replacing (62) into (60), one obtains 

O" + 3<H ( 1 + cfj O' - c 2 s AO + \m' + ( 1 + 3c 2 s ) ("K 2 - 7C)] O = 4nG N a 2 r5S , (63) 

which can be understood as the general relativistic version of Eq. (39). 

Finally, this evolution equation can be made to a much simpler, intuitive and tractable 
form: by setting 



Aa 2 3 ( H' 7C 

u = -— O with 6= J— =- and T^l -r + — r, (64) 

3 <H V2a 2 r «K 2 <H 2 

one can check after a few tedious but straightforward calculation that Eq. (63) takes the 
wavelike form 



- 2 2 Q"\ SnG N a 4 6 

+ C K -\U= — —TOS, 



3 <H 



(65) 



where I have replaced the Laplacian A —> —k in Fourier space. When there is no en- 
tropy perturbation (adiabatic perturbations), this equation is simply that of a parametric 



oscillator; the entropy contribution can then be seen as a forcing term. As it turns out to 
be the same as the gravitational wave case, I now move to those. 



THE EXAMPLE OF TENSOR MODES 

Tensor modes, being gauge-invariant from the outset, are free of all gauge-fixing sub- 
tleties, and can be computed straightforwardly. Einstein equations for those read 

E'l { + 2<HE' H + (2<K - A) E u = %nG n a 2 pn kh (66) 

where the anisotropic stress nu is usually set to zero, in agreement with the observations. 
Moreover, as we have seen, the spatial section curvature is also measured to be quite 
small, so we can safely set it to zero as well. Since inflation also set both these quantities 
to vanishingly (exponentially) small values, we have both observational and theoretical 
good reasons to set nti —> and 'K —> 0. 



Flat space polarization 

In order to understand what a tensor mode is, it is simpler to first consider the non 
expanding case in which we set 'H —* 0, so the Einstein equation for reduces to the 
wave equation 

aEij = (67) 

whose solutions I now discuss. 



Polarization. 

Let us consider for simplicity a mode propagating along the z direction, and pick the 
simplest possible solution of (67), i.e. oc cos[k(z-t)]. Now what is missing in this 
solution is the set of indices, which account for the polarizations. With ki = (0,0, k), 
the transverse and traceless conditions for En read E xz = E yz = E zz = 0, E xx = E yy and 



E xy = E yx . We are thus left with two independent solution, E xx and Exy say. These are 

the functions behaving as sines and cosines. 

The full solution can be expressed in terms of these functions together with a set of 

polarization tensors s + . and e x ., namely 
l j l J 
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(68) 



whose names stem from their effect on a test particle. 
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FIGURE 6. Left: Geodesic deviation rf of two geodesies defined by ^(/l)^ and s and represent- 
ing the trajectories of two test particles. Right: When a gravitational wave mode e + or e x passes through 
a ring of such test particles, they evolve as shown, producing the '+' or 'x' shapes with time, for different 
values of time in units of the period T. 



Observing a gravitational wave 

Let us consider two such neighboring test particles following their own paths x^(A)\ 

and (see Fig. 6) and assume they are originally at rest. The connections F 00 

being vanishing at first order in perturbations (only Ejj is present), as the wave passes, a 
particle originally at rest remains apparently so: the particle is moving with the reference 
frame. However, the perturbed curvature is non vanishing, so the relative geodesic 
motion is affected by the wave. The geodesic deviation = dx^/ds between these 
geodesies feels an acceleration given by 



dV 



= u a V a (uPVpn ti ) = R 



41 R u v u a nP, 

yap ' 



(69) 



so that the distance between the ring-forming particles changes with time as 

leading to the time evolutions shown in Fig. 6. 

The very simple cosine and sine solutions are obtained in the flat Minkowski case, 
and can readily be generalized to the expanding case with a scale factor increasing as a 
power law: they are then replaced by Bessel functions, see below. 



Cosmological gravitational waves 



Even under the simplifying assumptions made below Eq. (66), one still needs take 
into account the tensorial nature of the modes. Because this equation is linear, it can 
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FIGURE 7. Notations for the gravitational wave propagation: the wave propagates along the direction 
of k, and the dyad e a (a = 1,2), in the plane orthogonal to k, is used to define the polarization tensor. 



easily be decomposed into Fourier modes in the form 

£ y = -7^372 S f d 3 kPf j (k)iu A (n,kW kx , (70) 
a (2tt) 3/2 fr( J 1 

where A is a polarization index (running from 1 to 2 because of the two degrees of 
freedom in the tensor modes), and the polarization tensor P A . can be given explicitly as 
follows. Note the factor I /a which has been put here for further convenience. 

Polarization modes 

Figure 7 displays the configuration for the wave vector k and the orthogonal plane in 
which one defines a dyad e a (a = 1,2) satisfying e a ■ ej 7 = 6 a b and k ■ e a = 0. Recall that the 
tensor mode Eij is transverse and traceless, which translates into the polarization tensor 
as the requirement 

fc ! >J = 0, and Pfj6 ij = 0, (71) 
(we have set 7C — > and thus can identify y !; — > 8 l i). One can check that the choice 

P { - } ) = {e l ) i {e l ) J -{e 2 ) i {e 2 ) J and Pf) = (e x ) i (e 2 )j + (e 2 ) i (e l ) j (72) 

satisfies all the constraints (71). 

Plugging the form (72) into Eq. (66) then shows that both quantities n\ satisfy the 
same differential equation. We then simply set ^=1,2 = (the index T standing for 



tensor, we will later have a similar variable with an index 'S' for the scalar case) and 
obtain 



(73) 



which is the prototypical wavelike equation obtained in cosmological perturbation the- 
ory. It is interesting to realize that it can be obtained by varying the Einstein-Hilbert 
action expanded to second order in perturbation, namely (for the general case including 
curvature) 



d 3 jcd?7 y/y 



(74) 



in which one recognizes the action of a time varying mass scalar field. This observation 
lies at the heart of the idea of setting quantum initial conditions, as I will explain later. 



Time development of a mode 

Let us go back to the original equation (66) in the case of an expanding Universe 
dominated by a perfect fluid with constant equation of state. In this case, we have seen 
that the scale factor behaves as a oc if for some value of v, and the Hubble function then 
takes the simple form *H = v/rj. Eq. (66) thus becomes 

dEn IvdEij 

-{ + -— 1 +E lJ = 0, (75) 
dx z x ax 

where x = krj. As announced earlier, this is a Bessel equation whose solutions have been 
studied in details since the beginning of the 19 th century. They are shown on Fig. 8 and 
read 

i 

(76) 



I 

Eij = x2- 



AijJ v ^{x) + BijN v _i_{x) 



where the tensors A ;J and 5 ;J are, of course, transverse and traceless. 

The solutions show two extreme regimes, called sub- and super-Hubble. They refer to 
the characteristic ratio of the wavelength (k) to the Hubble scale (oc 77" 1 for a power-law 
scale factor), i.e. the variable x = krj. Long wavelengths (x «: 1) are strongly damped 
by the expansion, and so any initial motion is rapidly erased by the expansion-induced 
friction and the mode behaves as a constant. For short wavelengths (x » 1) on the other 
hand, the expansion is negligible and the mode behaves essentially as in a Minkowski 
universe: we recover the oscillations obtained earlier. 



Hubble vs Horizon 



At this stage, I think it is important to make a short comment on a commonly used and 
very misleading phrase, namely the use of "sub-horizon" and "super-horizon" modes 




FIGURE 8. Solutions of the gravitational wave modes equation (75) for a given polarization mode 
Ejj as a function of x = krj for radiation (v = 1, left panel) and matter (v = 2, right panel) domination. The 
dotted line represent the Hubble crossing time £77 = 1 : for £77 <K 1, i.e. deep in the super-Hubble regime, the 
solution behaves as a constant, while for the sub-Hubble case krj » 1, the solution exhibits the oscillatory 
behavior already encountered in the Minkowski case. 



instead of sub-Hubble and super-Hubble. As shown by Eq. (31), the horizon is a global 
quantity which depends on the entire history of the Universe. Early models were based 
on the singular big -bang followed by radiation and matter dominations, and it is easy 
to convince oneself that in this case, the integrated quantity d H is, up to an irrelevant 
numerical factor, given by the Hubble scale. In such a context, comparing wavelengths 
with the Hubble radius would indeed be equivalent to comparing them to the size of 
the horizon... but it is then not very clear what the meaning of these would have been! 
Indeed, in GR, modes larger than the horizon are actually not well-defined in the sense 
that setting initial conditions for them would be a direct violation of causality. 

In fact, one often reads that the modes are "frozen" because of some "causality" 
reason, with the meaning that a mode larger than the horizon could not evolve at all 
because of causality, as both ends of the mode would need to propagate faster than 
light to communicate, which is forbidden. I do not know what is the meaning of such 
arguments, and strongly suspect they have none whatsoever. A given mode consists of 
a linear combination of the two independent functions solving a second order linear 
equation 4 , with coefficients provided by the initial conditions. Then, for a super-Hubble 
wavelength, what happens is that the expansion very rapidly suppresses one of the 



This discussion is also valid in the case of scalar perturbations, so I do not restrict attention here to the 
tensor case. 



solution relative to the other, and one is left with the constant mode as discussed above. 
But this is in no way related with causality, on the contrary, it is a purely dynamical 
statement. 

We have discussed this point in more details in Ref. [7] for the specific case of scalar 
perturbations in bouncing models. 



Radiation-to-matter transition 

With the mode evolution known in any given epoch and a primordial spectrum, 
one should in principle be able to predict the observed spectrum. As it turns out, the 
theories that agree with the data predict an almost scale-invariant initial spectrum, i.e. 
one in which no particular scale is singled out. On the other hand, we know that such a 
scale should be present somehow, because the Universe, which was at very early times 
dominated by radiation, transitioned to the matter era 5 . 

The transition can be treated simply by introducing a new variable y = a/a eq , where 
a eq is the value of the scale factor at equality between radiation and matter, shown in 

Fig. 4 and defined through p m («eq) = Pr(«eq)- Given that p m = p^a~ 3 and p r = p®a~ 4 , 

we find that a eq = p^/Pr ' and finally that y = p m /p T . I leave as an exercise to the reader 
to show that the total equation of state w, defined as the ratio of the total pressure by the 

total energy density, is w = j (1 + y) _1 and the sound velocity is c 1 = | (l + |yj . 
The Friedmann equation (21) takes the form 

n 2 = -^pa 2 = -^p t (l + y)y 2 al q => <H 2 q = —^?L => * = 

(77) 

thus allowing to switch to the variable y whenever one encounters < H. Noting that the 
derivatives with respect to rj and y satisfy 

d dy d a' d ^ d 
drj drj dy a eq dy dy ' 

and defining the wavenumber characteristic of equality as k eq = < H< tq = a eq H eq , we find 
that Eq. (66) takes the following form 

^U^^ + ^iA.Eu^ (78) 



dy 2 2y(l+y) dy \k eq ) l+y 

whose analytic solution is not known... but we can solve it numerically for different 
values of k. This is done in Fig. 9. 



5 There was another transition more recently when the Universe became dominated by the cosmological 
constant or whatever it is which mimics it nowadays; I will not discuss this any further, but in principle, it 
could well lead to another scale in the data indeed. 




Iny 



FIGURE 9. Solutions of the gravitational wave modes equation (78) during the radiation-to-matter 
transition for 3 different values of the wavenumber k as function of the rescaled scale factor variable 
y = a/a eq . Long wavelength modes (small k) see the transition later than short wavelength modes (large k): 
the latter are more damped than the former, and the characteristic scale of the transition is thus imprinted 
into the subsequent spectrum. 



Let me conclude this section by returning to the form (73) of the mode equation. Its 
formal solution is known and given by the Born expansion, namely 



H T (k,r}) = a(ri) 



n dr , n da r a , 

A l (k)+A 2 (k) ^—+k 2 -5— dra 2 (T)fi T (k,r) 
J a L (j) J a z (cr) J 



(79) 



i.e., we have replaced a differential equation by an integral one! However, we have 
gained something in the process because the solution is known exactly in the long and 
short wavelength limits. They are, for short scales 



r » 



H T = A(k)e ih > + B(k)e 



-ikq 



(80) 



and for the large scales 



r « — 
a 



/u T = C(k)a + D(k) 



H dr 
a i « 2 (r)' 



(81) 



so that it suffices to plug one of these solution into (79) to obtain an expansion in either 
large or short wavenumber s. 

One final point, regarding the action (74) and the solution (80) valid for small wave- 
lengths, i.e. when the expansion can be discarded: in this case, we can consider fi T as 
an actual scalar field in a Minkowski universe, a system which we know how to quan- 
tize. Once the field itself is quantized, one can impose a given physical quantum state, 



for instance the vacuum. But this actually fixes unambiguously the coefficients A(k) and 
B(k)... and then, the theory becomes predictive! I shall return to this in the last part of 
these lectures. 

After having discussed the tensor modes in quite some depth, let me now sketch the 
scalar case which, although currently the only one observed, is also sufficiently more 
involved to require a special treatment of its own. This is done in Refs. [2] and [5]. 



DENSITY FLUCTUATIONS AND THE POWER SPECTRUM 

Since the tensor modes have not been observed yet, let us move on to the scalar case. 
For this, I will set a much oversimplified (and already disproved by the data!) model 
having 'K = 0, i.e. £jQ a = Q to t = 1, and assume all there is to consist of matter now (i.e. 
= 1) w i tn a currently negligible amount of radiation fl{? «; 1. 



Basic equations 



Before we move on to our specific example involving basically only matter (w m ) and 
radiation (w T ), let us write down the more general set of equations for a fluid having ./V 
constituents, e.g. radiation, dust, neutrinos, dark energy, dark matter, and whatever else 
a theoretician's brain can come up with. 

The total energy density p = J^ a p a and pressure p = YiaPa allow to define a global 
velocity through (p + p)v l = Yia^Pa +Pa)v l a . The total equation of state w = p/p and 
sound velocity cl = p'/p' can be obtained as 



Qw = 2^ &aWa and Q.c s = ^ 1 + & a c a 



(82) 



where each fluid sound speed is c 2 a = p' a /p' a . 

Similar calculations can be made at the perturbation level, yielding CIS = Yua Qa8 a for 
the density fluctuations and (1 +w)Q.v = Jj a ( \ +Wa)0. a v a for the velocities. The total 
entropy perturbation can be derived in much the same way as for the single fluid case, 
namely recalling that we set tSS = PT, we get wT = (dp - c^Sp} /p, and finally 



(83) 



showing that even if each individual fluid has vanishing self entropy perturbation (i.e. 
even if all Y a — » 0), the total fluid mixing entropy can be non vanishing. 

In principle, if the fluids are coupled, one should not necessarily assume them to be 
each independently conserved, but rather to satisfy (the condition in parenthesis being 
redundant) 



v A' i a — \la 



with 



and therefore 



7T = 



(84) 



assuming some action/reaction principle for the various fluid components. For the back- 
ground, setting <2^ = (-aQ a ,0), we have the generalization of (20) to a many-component 
fluid, namely 

p' a + 3<H(l+ Wa )p a =aQ a with J]Q a = 0. (85) 

a 

In practice however, since we shall here restrict attention to matter and radiation, we set 
the forces acting on the fluids Q v a —> 0. 

In terms of these variables, the relevant Einstein equation reads 

AO = ^«K 2 ^^X, (86) 

a 

showing how to relate the large-scale structure distribution (the density fluctuations) to 
the metric perturbations. The perturbed densities and velocities, when both the forces 
Q a —> and the self entropies F a —> are vanishing, follow the continuity and Euler 
equations 

(t^— ) +AV fl -3O'=0 and V' a + <HV a + ®+ -^—6% = 0. (87) 

Let us specialize for now on to the case of two fluids. Introducing the gauge-invariant 
relative velocity v and entropy perturbations S 

v = v a -v b and S = — , (88) 

1 + w a l+w h 

relation which can be inverted through 

+ ^^W = - 1 +n b s, (89) 



l + W a l+WbJ l + Wb 

the continuity equation can easily be restated as 



S' = -Av-3<Hf, where f^^^-^^, (90) 

1 + W a l+Wb 



while Euler equation reads 

v' = -m-(cl-cl)-^ + 

v °> 1 +w 



C a (y+W b ) — +C Z b {\+Wa) — 



l+W 



-Tab. (91) 



The basic idea now consists in solving Eqs. (90) and (91) together with (63) so as to get 
a complete solution for the distribution of 6 C now through (86): this density perturbation 
spectrum can then be directly observed as the large-scale structure distribution. This is 
more easily said than done, and to begin with, one needs to impose initial conditions, to 
which I now turn. 



Adiabatic and isocurvature initial conditions 



In order to impose a complete set of initial conditions, we need to know the number of 
independent degrees of freedom. As soon as one knows all the fluid variables, the system 
is fixed, namely knowledge of all the 6 a and v a is enough. For the two constituent fluid 
we are dealing with, this means we need 4 independent conditions. With the previous 
variables, we can re-express all these in terms of the sums 6 C and V, so that then Eqs. (86) 
and (61) provide O and O' respectively. We are then left with the relative velocity and 
entropy perturbations v and S . 

There are basically two sets of initial conditions which are used, the so-called adia- 
batic and isocurvatures ones. They are defined by the following conditions. 

• Adiabatic initial conditions: the entropy perturbation (88) vanishes at the initial 
time, while the Bardeen gravitational potential O is a constant, so we have 

S = => = and O' = 0, (92) 

1 + W a l+W b 

leaving two arbitrary initial numbers, <D and S' say, to decide of the forthcoming 
mode evolution. 

• Isocurvature initial conditions: the opposite, and complementary, situation consists 
in setting initial conditions such that there is no initial metric perturbation, i.e. we 
set at the initial time 

O = and 6 C = => ^Q fl c5£=0, (93) 

a 

which, again, leaves 2 arbitrary numbers to be set, for instance the values of <E>' and 
the initial entropy S . 

These conditions essentially reflects all the possibilities, and the "real" initial pertur- 
bation should be a linear superposition of those. 

It should be reminded at this stage that these initial conditions ought to be set at the 
point in time after which we can evolve them with sufficiently precise knowledge of 
the cosmic history. Normally, this is done using a simulation code taking into account 
all known relevant cosmological phenomena. This means we assume that the initial 
spectrum of perturbations is propagated through the almost entire history of the Universe 
from this original time... the question then remains of what is this initial time, and how 
can we even suppose we know anything at all about it? That will be the subject of the 
final section, but for now on, let us concentrate to the actual evolution and the spectrum 
we might get now so as to be able to compare with observational data! 



Mode history 



As discussed above, we now specifically restrict attention to a flat, matter-dominated 
Universe containing a tiny amount of radiation, so that Eqs. (87) now read, in Fourier 



space, 



<=fc 2 V m + 30' and V^ + W m + O = 0, 

4 1 

Sf = -k 2 V r + 40' and V' m + O + -d? = 0, 



(94) 
(95) 



and we assume we know, somehow, the initial conditions for the perturbations deep into 
the radiation epoch. Adding the Fourier-transformed Poisson equation 



-k 2 ® = -<H 2 
2 



Q m «5»+Q r ^-3<H 
closes the system which we now solve. 



(v m + ^v r 



(96) 



Initial conditions in the early radiation epoch 

The Universe will have to go through the radiation-to-matter transition, and we thus 
switch to the relevant time variable y = a/a eq defined above Eq. (77) and in terms of 
which the relative density parameters read 

y 1 
£l m = and Q, = 



l+y 1 l+y 

Since now S = 6 m - |(5 r and noting the relationship between the comoving total density 
perturbation and the Bardeen potential 

, 4 / k \ 2 y 2 
6 C = -- — O, 



3\k sq ) l+y 

we see that the entire system reduces to the set 

' d 2 1 / 1 8 \ dO O _ 2 / c yS 

dy^~ + 2y\ ~T+y~ + 4 + 3y)~dy + y(l+y)(4 + 3y) ~ y2(4 + 3y)\ ~T+yJ' 

d 2 S 3y + 2 dS Ilk \ 2 / yS 



(97) 



{ dy 2 2y(l+y)dy 4 + 3y\fc eq / \ l+y 



Once we have the solution to the system (97), we can reconstruct the density perturba- 
tions through 

d<D\ l(l+y)6 + S (l+y)S-yS 

^=^-2\d> + y—\ and 6 m = ± , 5 r = K y) y . (98) 

dy) l + |y l + y 



We now want to impose initial condition very early on, when the Universe is radia- 
tion dominated, i.e. for yi n i «: 1, and we are interested in cosmologically relevant wave- 
lengths, i.e. those satisfying k «: r H 1 m. If we decide for adiabatic initial conditions, this 



means we can demand 



dO 

= ini , and — 
dy 



turning Eq. (98) into 



y=yh 



Si ~dy 



y=yim 



0, 



(99) 



1111 3 \<H m 



leading to 



Oiri, £ c . =6 C . and £ c . . = -6 C ., 

1111 ' r,im mi m,ini a inr 



(100) 



cJ N . . =£ N . and £ N . . = -£ N .. (101) 

r,ini mi m,im A mi v ' 



These initial conditions mean that the density ratios are constant everywhere on the 
initial hypersurface: both density perturbations behave in the same way everywhere, as 
illustrated in the left panel of Fig. 10. 



Space point x 




Space point 



FIGURE 10. Adiabatic and isocurvature initial conditions in terms of the primordial density fluctua- 
tions: the left panel shows the adiabatic condition, where the radiation (full line) and matter (dashed line) 
density perturbations are everywhere following the same pattern of fluctuations, while the isocurvature 
condition, represented on the right panel, has a dominant contribution coming from the matter fluctuation 
(full line) together with a negligible amount of radiation density perturbation (dashed line). 

Setting isocurvature initial conditions on the other hand amounts, in this case, to 
imposing at the initial time the relations 



S =Si 



and 



dS_ 
dy~ 



dO 

= Oi = — 

y=yim ay 



y=Tm 



= 0, 



(102) 



which translates, in terms of density perturbations, into 



6 C . = 0, 
mi ' 



rC 



r,ini J 1111 m,ini 1111 



= Si 



and S C .=6 N . 

1111 mi 



(103) 



with all velocities vanishing. This is also illustrated in Fig. 10, on the right panel. 



Transfer function 



Let me now move on to the transfer function, which is defined as the ratio of the 
observed spectrum of perturbations now, i.e. essentially the large scale structure distri- 
bution in the sky, with the primordial spectrum, as calculated by high energy physics; 
this primordial spectrum is the topic of the final section. 

Let us consider a given mode, i.e. a given wavelength k, and solve the equations of 
evolution (97). Without entering unnecessary details that can be found elsewhere [2, 5], 
suffice it to say that the main behavior of the perturbations depend on two quantities, 
namely the time at which they are evaluated relative to the equality time rj eq , and the 
comoving wavenumber of the perturbation, again relative to the value k eq , defined above 
Eq. (78). Fig. 1 1 shows the various possibilities by comparing the wavelength with the 
Hubble factor entering the evolution equations. 



FIGURE 11. Evolution regimes for perturbation modes. All modes initiate at time T]i a i during the 
radiation dominated epoch, and are all super Hubble initially, i.e. they all satisfy k < Him- Depending 
on whether the comoving wavelength is larger or smaller than the Hubble scale at equality [recall 
&eq = 'Heq = ^(^eq)], tne mode becomes sub Hubble before or after the equality rj e0i . This leads to different 
time evolutions, and a final spectrum that, more or less independently of the initial spectrum, will have 
the equality scale k &q imprinted in it. 

One finds the following time developments: 




Radiation 
Domination 



Matter 
Domination 



• k<k 



eq 



• k>k { 




- 7] m i <7]<l/t.®~ O ini -» 6 m OC T] 2 , 



- l/k<r]<r] eq : 6 m oc\na, 

- 7] eq < TJ < T]Q\ 6 m OC J] 2 . 



We see that for most of the time, the density perturbation in the matter fluid evolves 
as the square of the conformal time, except for the modes whose wavelength is smaller 
than the equality scale: those become sub Hubble during the radiation dominated phase, 
during which they cannot grow because of the photon pressure. This is shown in Fig. 12. 



FIGURE 12. Time evolution of a mode with the corresponding transfer function as a function of 
the comoving wavenumber k. The density fluctuation for modes becoming sub Hubble during matter 
domination and having k < k eq essentially grow at all times as oc rj 1 , corresponding in fact to different 
behaviors with the scale factor. For a finite amount of time however, modes becoming sub Hubble during 
radiation domination cannot grow before the advent of matter domination and thus acquire a logarithmic, 
i.e. almost constant, time behavior. As a result, their amplitude increases less than the other mode's 
amplitude, leading to a transfer function T(k) as indicated in the insert. 

The previous time evolution transforms into a change in the spectrum for scales above 
k eq . Indeed, a mode 6 K with k < k eq evolves essentially as rj 2 all along, so that 6q ~ S^ tj^, 
with the index "0" still denoting the present-day time. Similarly, a mode 6 > with k > k eq 
evolves as rj 2 only up until 77 ~ k , at which point it behaves roughly as a constant. As a 
result, when it starts growing again, at T] sq , its value is 6 > (rj sq ) ~ S^-k iji. The transfer 
function T(k) is now defined as the ratio between the final (evaluated at 770) an d initial 
(at 77i n i) density perturbations, namely 



and the calculation above shows that T(k) ~ 1 for k < k eq (long wavelengths) and 
T(k) ~ k~ 2 for k > k eq (short wavelengths). The insert of Fig. 12 also shows the typical 
behavior of the transfert function. 
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6(k,r]o) = T(k)6(k,T] mi ), 



(104) 



Perturbation spectrum 



With the transfer function known, we can now derive the actual large scale structure 
distribution, can we? Well, in fact not quite yet, as there is something missing: the initial 
distribution 6(k,rji n i). This will be the subject of the last section, as I already mentioned 
a few times, but something can already be said at this stage, in particular by looking at 
the data. 

The only analysis that can be done of all the observation is statistical in nature, as 
we now understand the actual density distribution to be but a particular realization of a 
statistical ensemble, so that the density field itself is now seen as a random variable at 
each point. What is actually measured then is the correlation function £(r) of the density 
field, defined by 

Hr) = (6(x)6(x + r)), (105) 

where the mean value should represent an ensemble average. In practice however, we 
have only access to one such realization, and we replace the ensemble average by 
a spatial average. Note at this stage that the cosmological principle implies that the 
distributions should be isotropic and homogeneous. As a result, the correlation function 
should only depend on the distance scale r and neither on the particular direction choice 
r/r nor on the specific point x. 

Moving to the Fourier space, we can write 

6{r) = f 7^^*>^ r - 0° 6 ) 
J (2tt) 3/2 

whose spectrum P$ (k) stems from the two-point function in Fourier space, namely 

(6(k)6(p))=Ps(k)6(k + p). (107) 

I leave as an exercise to show that it is the Fourier transform of the correlation function, 
i.e. 

<P ^) = f -^£(r)e* r , (108) 

and, as expected again from the cosmological principle, it also does not depend on the 
direction k/k but merely on the wavenumber k. 

In a way, the power spectrum can roughly be seen as the square of the density 
distribution. Therefore, we also have the relation 

P 6 (k,r]o) = T 2 (k)Ps(k,i]M), (109) 

which is the equivalent of (104) for the spectra. 

We shall see later that the expected primordial spectrum actually scales like 
P6(k,T]ini) °c k, and so the observed distribution should scale as k for long wave- 
lengths where the transfer function is independent of scale, and as k~ 3 for shorter 
wavelengths. Fig. 13 roughly confirms these expectations. 

A quick estimate of the equality scale k eq is provided by 

k eq = H a ^2Q? m (l + Zeq ) ~ 0.072Q> 2 Mp C - 1 A eq = k£ ~ 
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FIGURE 13. Actual observations of the large scale structure distribution showing a linear behavior in 
the wavenumber for large scales, followed by a decrease as £~ 3 for small scales. The spectrum peaks 
around a few hundreds h~ l Mpc, which thus indicates the value of the equality scale, i.e. the Hubble 
radius at j] eq . 

which is estimated to a hundred Mpc, falling a bit short of the actual value. This is due 
to our very rough approximation according to which our universe model only contains 
matter and radiation. 

The actual transfer function is in fact much more complicated to calculate once one 
takes into account all the cosmological constituents. For instance, if there is a so-called 
hot dark matter component, consisting in relativistic degrees of freedom at the time of 
decoupling, e.g. neutrinos, it has the effect of streaming easily throughout any initial 
perturbation, thus wiping out very large scales. These scales cannot grow anymore, and 
this produces an exponential cutoff in the transfer function. 

INITIAL CONDITIONS: QUANTUM VACUUM FLUCTUATIONS 

So far, this presentation contained essentially no new physics, merely basic applications 
of general relativity and fluid dynamics. It can all be made much more precise, e.g. to 
include thermodynamics and using the Boltzmann equation, but this would also be well- 
known physics and in no way can ever provide what we are seeking, namely the initial 
conditions for the perturbations we have just calculated. In fact, measuring the large 



scale structure distribution while knowing all what precedes is akin to measuring the 
initial conditions, which is kind of useless if we don't have a theory that predicts them. 

It turns out that soon after the advent of inflation, which was originally the first 
convincing scenario that was proposed to cure the puzzles discussed in the first section 
of these lectures, it was realized that the accelerating epoch had the ability not only 
to enhance already-existing perturbations, but also to produce those when quantum 
vacuum fluctuations were taken into account. Since that time, other models, for instance 
including an initial contracting phase and a bounce to connect to our currently expanding 
epoch, have been devised that also provide the required initial conditions, and they can 
be tested quite accurately with the more and more precise data that are accumulating. 

Although the inflationary paradigm is the best accepted one to describe the primordial 
epoch, I would like first to emphasize that it is not established beyond any doubt (as one 
sometimes reads!), so that looking for challengers is still a reasonable activity. However, 
in what follows, I will restrict attention to the inflationary case as it is easier to implement 
and pedagogically more convenient. 



Back to the background 



Inflation provides explanations to the standard model puzzles by means of accelera- 
tion of the scale factor, namely for a finite but sufficiently long period of time, we have 
a > 0. Eq. (18) then implies that the pressure should be more negative than a third of the 
energy density (which is always assumed positive). This is easily achieved by means of 
a slowly-rolling scalar field ip whose action we take to be 



with a yet- undefined potential V((p) 



1 



(d<py + v(<p) 



-gd 4 x, 



(110) 



Slow-roll parameters 
The stress-energy tensor derivable from the action (1 10) is 



1 



■{dipY + Viy) 



1 



gfiv 



p=-<p 2 + V and p = -<p z -V, (111) 



■ 2 



where the definition of the energy density and pressure assume the field <p to depend only 
on time in order to satisfy the background symmetries. We see that the r.h.s. of Eq. (18) 
reads p + 3p = 2{p 2 - V), which can be negative quite easily provided the kinetic term (p 1 
is sufficiently small compared with the potential. Because then the velocity of the field 
is tiny, this is why we speak of slow-roll phase. 



The Einstein and Klein-Gordon equations then transform into 



8ttG. 



3 

SnG 



H 1 

a 
a 

ip + 3H<p + — 
cap 



'K 

a 2 ' 



0, 



(112) 



where the last of these, merely reflecting the conservation of (1 1 1), is not independent of 
the first two. Combining those actually yields H = -4nG N (p 2 + *K /a 2 : as (p 2 is assumed 
small, the natural tendency for H is to decrease as the scale factor increases. But this 
makes the scale factor increase even more rapidly, so the spatial curvature term becomes 
more and more negligible. In fact, this is an attractor of this system of equations, and 
therefore, for now on, we will assume 'K —* with the meaning that spatial curvature 
terms are exponentially smaller than any other. 

Applying the slow-roll conditions Cp 2 «; V and ip «; 3H(p, we find the relations 
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which are consistent with the original assumptions only if \H\/H 2 «; 3/2. More gener- 
ally, one can define two small parameters, called the slow-roll parameters, by 
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Inflation goes on for as long as a > 0, which translates into s < 1 . 

The simplest solution for this model consists in demanding s to be constant; in this 
case, the scale factor can be calculated as follows. First, I recall the relationship between 
conformal and cosmic time, namely 



dt = adr; => 

so that, noting we also have 

,M\_J_da ldH 
\aH)~ Ha 2+ aH 2 ~ ' 

this means that the relation 
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should hold, and if s is roughly constant, we can perform the integration, leading to 



which I can invert and obtain that the scale factor behaves as 



i.e. we have an exponential quasi de Sitter phase. The parameters to and t]q are con- 
stant of integration necessary to pass from conformal to cosmic time. Eq. (117) shows 
moreover that inflation occurs in the regime where rj < and 77 —> 0~ . 



Two explicit examples 

The simplest example one can think of is that for which the scalar field is merely a 
massive free (non interacting) field, namely the potential reads V((p) = jm 2 (p 2 . In this 
case, the system (113) reads 

mM m 
3H(p + m 2 (p = => < p(t) = <pi Bi -——^t, (118) 



where use has been made of 



The scale factor is then 
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a (t) = fl ini exp I ^ [p 2 m - ip 1 (t)] J , 

from which one obtain the slow-roll parameters as 

M 2 

s=—^ and 6 = 0. (119) 
Amp 1 

As e varies with time, we can easily calculate when the slow-roll phase ends, namely 

for (p = (ff = M pl / V47T, so the number of e-folds of inflation is N = 27T^(p m [/M p ^j - \ . 
In order to solve the cosmological puzzles, we know that we must impose /V > 70, 
leading to the requirement that the initial value of the scalar field should be of order 
<£i n i - 3M P1 . One might think that this could be a problem, as such a high energy scale 
would require quantum gravity to be described properly, and of course we do not have 
such a theory. However, what actually matters is not the field value itself, but the energy 
density that it stores. Under the slow-roll hypothesis, this means the potential energy, 

i.e. Vim ~ I (mM n ) «: provided the scalar field mass is much less than the Planck 
scale. As we shall see below, this is exactly what is required from the data. 

Another useful model is the so-called power-law inflation, for which one demands 
the scale factor to increase as a power-law instead of an exponential, while still being 
accelerated. Explicitly, this is 

a = a„(-77) 1+/J a = a t t p with 1 +8= (120) 

l-p 



where a n and a t are constants, and p > 1 to ensure that a > 0. 

Integrating (113) equations again, one obtains the scalar field behavior 
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and the potential it evolves in 



V(cp) = Viniexp 



together with the slow-roll parameters: s = 6 = I /p. Since there is no time evolution in 
this case, we see that such a model has merely a pedagogical value, as inflation never 
ends in this case. However, it can really be useful because many features are calculable 
in an analytic way. 

Having settled and somehow implemented the inflationary phase, let us see what 
happens to fields living in such a background. 



A test field in inflationary background 

We shall here follow the evolution of the simplest case, namely that of a test field in an 
inflationary background which, to make things even simpler, we shall assume takes the 
form of a quasi de Sitter expansion, i.e. an actually exponential expansion, or Eq. (1 17). 



Massless scalar field 

Let us begin with yet another simplifying assumption, namely that the test scalar field 
X is massless, so its potential vanishes. The Klein-Gordon in the expanding background 
reads 

X + 3Hx-\&x = => x" + 2W + k 2 X = 0, (121) 

a 1 

where in the last equality I switched to the conformal time and took the Fourier trans- 
form. Now, setting v = ax, thereby defining v, we find that 

v" + [& 2 -(<H'+<H 2 )]v = Qr 1 // + | Jfe 2_^j v = l // + | Jk 2_^fj v = . (122) 

where I have used Eq. (1 17). 

Eq. (122) is the same as that driving the evolution of the tensor modes (73), and is of 
the generic form of a "time-independent" Schrodinger equation in a potential (regarding 
the conformal time variable as the equivalent of the spatial coordinate): 

h 2 d 2 iA d 2 v r ? i 

- — -^ + [V(x)-E]^ = <=> — + \k 2 -U(r ] )]v = 0, (123) 
2m dx z dr] z L J 



provided one identifies V(x) with U(rj), E with k 2 and rescale everything to cancel out 
the -h 2 /(2m). In fact, whenever there is no entropy perturbation, it is the most generic 
form we will ever encounter. So for now on, I will assume a generic potential U (q), and 
discuss the actual mode evolution. 



Evolution regimes 



The potential during inflation grows like rj~ 2 when rj — » (T, but this is merely an 
artifact of our approximation. In a realistic scenario however, the potential might look 
like that represented in Fig. 14: it starts growing during the phase of inflation (or 
in general any such phase during which primordial perturbations are produced), and 
reaches a maximum, after which it decays. These last phases would usually represent 
radiation or matter domination, at which times we would observe the mode somehow: 
the scale factor would then behave as a oc rf, so the typical potential should look like 
a ,, /a=j3(fi-l)r } - 2 . 



Regime of quantum fluctuation amplification 
(Inflation, bounce...) 



Radiation domination 
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FIGURE 14. Different regimes for the time evolution of a scalar mode: the potential in Eq. (123) 
starts growing initially during the perturbation production epoch, then stabilizes for instance at the end of 
inflation or at a bouncing point, and then decays again while getting into the radiation or matter domination 
era (or any other relevant subsequent regime). For a short wavelength A sw , the wavenumber k 2 is at all 
time larger than the potential, which therefore doesn't affect the mode evolution: Eq. (123) then indicates 
a simple oscillating behavior at all times. On the other hand, for a larger wavelength, A LW , i.e. smaller 
k 2 , different regimes can be identified: initially, in region I, the mode oscillates as k 2 » U(rj), then there 
is a transition through region II in which the mode passes below the potential. Then in region III, one 
is in the opposite situation where k 2 <k U(rj), and the mode now consists in a growing and a decaying 
solution. Finally, region IV connects to the standard cosmology, the mode is above the potential again, 
and therefore starts oscillating again; these oscillations are those one observes in the Cosmic Microwave 
Background which I did not have space to discuss here. 



Figure 14 summarizes the discussion of the tensor modes, with the same kind of 
solutions (80) and (81); I shall not repeat this analysis here, but suffice it to say that it 
also applies to most known cases as very often the potential has the form of the second 
time derivative of a function over this function. 

In the special case of de Sitter expansion, i.e. (1 17) with e —* 0, the solution is known, 
since this is then a quite simple Bessel equation, and we have 



Vk 



(r ] )=A(k)e-' kr >ll + J-) + fl(*)e*»(l - -L), (124) 
\ ikrj) \ ikrjj 



where A(k) and B(k) are yet-unknown function depending only on the scale k. 

The massive scalar field case can be obtained in a very similar way as it suffices to 
replace (123) by 

r2 m 2 /H 2 -2^ 



whose solution is again another linear superposition of Bessel functions of index v, with 

9 

4 



v 2 = 9 m 2 /H 2 



Quantization 

All what precedes does still not tell us what initial conditions we should use, or, in 
other words, given (124), what should we take as functions A(k) and B{k)l 

To achieve this goal, we need to quantize our system, which is quite simply done 
when we have discussed the action expanded to second order. 



Expanding the action 
The action for our scalar field, still without a potential to keep things simple, is 

S = j \ {d X f V=^ 4 * = \ ja 4 \- X ' 2 + (Vx) 2 } d\ (125) 
which we can express in terms of the variable v as 
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variable mass scalar field in Minkowski space surface term, irrelevant 

showing it is nothing but the action for a simple scalar field in Minkowski space... usual 
technique of quantum field theory can now be applied, and we will have the possibility 
of choosing a specific quantum state to provide the initial conditions. 



Canonical quantization 
We can expand the field v as any standard quantum field through 

[v k (rjW kx a k + v* (i7)e-**4] , (126) 

the second term being the hermitian conjugate of the first. 

Quantization is achieved by promoting —* to an operator in the Fock space of 
field configurations and imposing the canonical commutation relations 

[a k ,a\}=£P\k-q). (127) 

These relations can be seen as stemming from the actual field quantization: defining the 
conjugate momentum 

n = — = v — > operator n, 
ov' 

with the Lagrangian being the integrand in the definition of the action, the Hamiltonian 
follows 

r 1 r I „tt\ 

44, 



H = J (v 'n - £) = X - J (n 2 + d.vd'v - ^ j d z 



and we can impose the standard equal time commutation relations for the field operators, 
namely 

[r(jc.;/),r(j.//)j=() = [ir(jc.//).ir(j,//)j and [Hx,T]),ft(y,T})] = i6 (3) (x-y). (128) 

These commutation rules are consistent with those above (127) only provided the Wron- 
skian W(k) = v^v'* - v*v' k is normalized to W = i since one gets directly from the field 
expansion 

d 3 k 



[Hx,rj),n(y,rj)]= [ ^e*^>W(Jfc) 
J (2n) 



(2ny 

We are almost done, having merely to define the relevant state to assume as initial 
condition. 



The vacuum state 

In quantum field theory and therefore here as well, the vacuum state is that which is 
annihilated by all the so-called "creation" operators a^, namely 

a k \Q) = for all k, 

and all other states are obtained by repeated application of the operators a£ on |0>. 



In the limit \krj\ » 1, i.e. for large negative conformal times where we indeed want 
to impose our initial conditions, we are back to the usual massless scalar field in a 
Minkowski space time, and we know that the vacuum state must therefore satisfy 

Vk — > —J=, 

as indicated in any standard textbook on quantum field theory. Given the previously 
obtained solution, this leads to the so-called Bunch-Davies vacuum state 

^07) = ^(l + ^)e-^, (129) 

which now provides a closed form initial solution for our perturbation. It is with such 
initial condition that one finally gets the scale-invariant spectrum which one compares 
with the observational data (and it works!). 

The power spectrum is now obtained from the 2-point correlation function £ v (x - y) = 
(0\v(x,r])v(y,rj)\0), which gives 



|v fc | 2 — , (130) 



after integration over the angles and setting r = \x-y\. 

It turns out that for large scales, i.e. super-Hubble modes, one finds that the properties 
of a quantum field are the same as that of a classical stochastic field with gaussian 
statistics. In particular, this means we can replace the quantum averages by statistical 
ensemble averages. For the stochastic variables, we find a power spectrum that scales as 

,,«-£,,<»- 

in other words a scale-invariant spectrum. 



Realistic perturbations 

If one wants to take into account all the actual complication of what is really going 
on, one needs to consider all scalar, vector and tensor modes of the metric and treat 
them including all possible effects. Although this is a very complicated task, some 
situations allow to say something however. For instance, assuming the scalar field (p 
to drive the early history of the Universe, one finds that the gauge-invariant degree of 
freedom generated by the variable 

tyV;£ = -. (132) 
Ti a 

thus defining the so-called Mukhanov-Sasaki variable v, is enough to describe all the 
fluctuations in a single field inflation. Expanding the action to second order in perturba- 



tion and getting rid of the surface terms just as above, one arrives at 
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(133) 



thus showing what I previously said, that the typical equation of motion is always of 
the same form. As the same analysis applies, one thus obtains a way to set up initial 
conditions by assuming quantum vacuum in the early stage of the Universe. This leads 
to a natural way to obtaining a scale-invariant spectrum that fits extremely well all the 
known data. 



CONCLUSION: CONDITIONS FOR ALTERNATIVE SCENARIOS, 

THE BOUNCING MODEL 

The cosmological scenario, as discussed in the notes above, represents a major achieve- 
ment in physics, performed in less than 100 years! In that time, it has been established 
that the Universe itself could be treated, studied and understood as a regular physical 
system, despite the fact that it seems to contradict, by its very uniqueness, the usual as- 
sumptions of the scientific method. In fact, cosmology somehow extended the inductive 
method, replacing for instance repetition of experiments by repetition of measurements 
in different directions, in other words, replacing ensemble averages by ergodicity. 

To summarize, we now have a rather clear view, basically, of what happened dur- 
ing the last 13.7 billion years, with detailed calculations comparing amazingly well with 
observations. Although I did not discuss them all, but these observations range from con- 
sequences stemming directly from nuclear physics (nucleosynthesis), thermodynamics, 
fluids mechanics, gravitational phenomena, and, as sketched in the last section above, 
the relationship between gravity and quantum physics! That only a bunch of "unpleas- 
ant" features are present in the data with the overall picture being generally consistent is 
absolutely astounding and very often not given enough emphasis. 

Now our cosmological model, precisely because of its successes, can be scrutinized 
with exquisite attention to unveil any possible new mechanism the we would not have 
thought about. This is how detailed examination of specific objects (Type la Super- 
Novae) and their redshift distribution revealed that the Universe appears to be currently 
accelerating (see however D. Wiltshire's contribution in this volume for an alternative 
understanding of the data), leading to a new component, dubbed dark energy, among the 
various fluids pervading the Universe. When added to the other components, it permits 
to fit all available data, including large scale structure distribution, SuperNovae or the 
Cosmic Microwave Background (CMB) fluctuations (see Fig. 15). 

There are many things lacking in this presentation, including not only the CMB fluc- 
tuations themselves, but also its polarization and measurement thereof, non gaussiani- 
ties, baryonic acoustic oscillations, and many others. One special point I would like to 
emphasize over however is that when all those data are taken into account, the whole 
thing can serve not only to check the currently accepted paradigms and models (an in- 
flationary phase followed by radiation, matter and cosmological constant dominations, 
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FIGURE 15. Comparison between various measurements of the CMB fluctuations, i.e. essentially the 
Fourier transform of the angular power spectrum of the microwave light coming from the recombination 
epoch (see other contributions in this volume). The figure shows the most recent nine-year WMAP data, 
together with those coming from SPT and ACT for the small scales. The standard ACDM model shown 
here is merely coming from a fit with only the WMAP data, which is then used to predict the higher 
multipole data: clearly a quite precise and correct prediction! 
Figure taken from [8] in which all references to other data are provided. 

the so-called ACDM model 6 ), but also to explore alternative possibilities: although the 
inflationary paradigms seems to provide satisfactory explanations to most of the cosmo- 
logical puzzles, it is not still proven beyond reasonable doubt, and besides, it has a few 
problems of its own. 

Inflation is based on well under control physics, i.e. GR and quantum/classical scalar 
fields, it can be implemented in high energy theories such as Grand Unification or string 
theories, and it makes predictions which have been experimentally shown to be compati- 
ble with observations... why would we therefore like to find any alternative at all? First of 
all, it does not really solve all the puzzles, as in particular, the question of homogeneity, 
although admittedly alleviated, is by no means solved. Moreover, an inflationary phase 
usually begins from a singularity, or from a quantum gravity fluctuating phase, which 
is not understood at all. Related to this is the fact that even the largest possible scales 
observed today, i.e. that of the Hubble radius, must have inflated and expanded from a 
time where it was actually smaller than the Planck length itself. Setting initial conditions 



" Meaning A for the cosmological constant, and "Cold Dark Matter", this phenomenological model 
describes with the minimal possible number of parameters the current data. 



there is, to say the least, debatable. Finally, providing challengers is always a very good 
way to test a theory, so inflation itself benefits from alternative models. 

Most alternative to inflation present, in one way or another, a contracting phase 
preceding the currently expanding one, to which it is related by means of a bounce. 
This is not a new idea, as it was in fact suggested in the 1930's by Tolman and Lemaitre, 
i.e. much before any inflationary scenario was even thought about. In course of time, 
bouncing scenarios were repeatedly proposed, as discussed in Ref. [9]. One might 
immediately argue that this seems to create more problems than it solves, since in 
particular it is very difficult to implement a bouncing phase in the framework of GR; 
however, the bouncing model also addresses different issues. For instance, there is 
of course no question of the primordial singularity, which is avoided by definition! 
Moreover, the horizon (31) can easily be made infinite if the initial time is sufficiently 
large and negative, i.e. in the limit t m { — » -<x>. Flatness is also quite a natural achievement 
of the bounce, as I discussed earlier. 

Now perturbation theory ought to be valid as well in a contracting background, so 
basically, all I said before applies straightforwardly in such a new framework. What 
needs be done then is to evolve similarly set vacuum initial conditions in the contracting 
Universe all the way to the bounce and up to now. In general, what happens is the 
following: whenever the relevant equation of the perturbations takes the form (123), the 
potential U(r]) can be more complicated than that shown in Fig. 14, and in particular 
it often happens that the term k 2 - U changes sign more than once or twice. As a 
result, the primordial spectrum starts oscillating before it gets amplified again, and one 
expects oscillations on top of the usual and expected oscillations. For the time being, no 
observation has been made along these lines, but one can hope to see those in the future, 
e.g. with Planck data. 

Finally, I should say that the perturbation question is somehow an open one in 
bouncing scenarios, and for many reasons. The first concerns for instance the vector 
modes: as I said before, one usually neglects them as they decay anyway with the 
expansion. Clearly, during contraction, one expects vector modes to grow, and therefore 
they might pile up to produce unwanted non linear vector-like objects, thus ruling out 
irremediably the corresponding model. Therefore, one needs to check every model and 
initial condition setup, although the situation is often quite unclear because without 
any specific coupling with the matter fields, the vector modes are not dynamical, so 
setting initial conditions for them is not feasible in any known natural way. Scalar modes 
themselves can grow very large, but then comes the question of gauge: is it absolutely 
clear that a large value for, say, the Bardeen potential, means the theory becomes 
non linear? As a matter of fact, this is yet undecided, and there are good arguments 
suggesting that providing there exists a set of variables that behave perturbatively all 
through the evolution of the Universe, then this set of variables should be used, at the 
cost of breaking "gauge invariance", and the theory would still make sense. 

Both inflationary and alternative models will probably be with us for still quite a 
while, unless some (always possible) unexpected prediction or observation comes in 
the way. In any case, we are living a very exiting period, not only of the history of the 
Universe itself, but also in cosmology where paradigm shifts are happening and new 
developments are proposed at an ever increasing rate. With the advent of forthcoming 
data (Planck of course, but also all the new proposals that just await actual construction), 



it should not take long before new ideas come in the front stage... hopefully, most of 
these notes will remain essentially valid. 
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